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We propose a new method to measure the coherence time of superconducting 
phase qubits based on the analysis of the magnetic-field dependent dc non- 
linear Andreev current across a high-resistance tunnel contact between the 
^ ', qubit and a dirty metal wire and derive a quantitative relation between the 

subgap I-V characteristic and the internal correlation function of the qubit. 
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a" 
o 

^ ■ 1. INTRODUCTION 



^ 



While new algorithms El'Bu are pushing the frontier of quantum compu- 
tation, their practical implementation requires novel ideas for the design of 
C^ \ qubits which can cope with the contradictory constraints of scalability and 

long decoherence time. Solid state electronics offers scalability but involves 
macroscopic elementary blocks which usually interact strongly with the en- 
vironment and, thus, suffer from a short decoherence time. In contrast, 
atomic scale designs (trapped atoms, photons in cavities, nuclear spins) can 
be easily decoupled from the environment but are difficult to combine into 
large and complex devices. The ideal system for a qubit implementation is 
a macroscopic device with two quantum states decoupled from the environ- 
ment which can be manufactured and combined with others by conventional 
techniques. This ideal is closely approximated by the superconductingDhase 
qubit (SCPQ) that consists of a frustrated Josephson junction circuitu with 
two states distinguished mostly by their superconducting phases. 

The simplest example of such a circuit is a low inductance SQUID loop 
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Andreev contact 



Fig. 1. SCPQ device with Andreev probe; the area Sq of the "intrinsic" loop 
of the qubit is equal to the area Sf of the NS-QUID loop. 



consisting of three (or more) Josephson junctions frustrated by a magnetic 
field that creates a flux <^ = $o/2 through the loopJ3 see Fig. 1. The two 
states of the qubit are defined by the phase changing clockwise or counter- 
clockwise around the loop. The advantage of an "all phase" design is that 
the superconducting phase is not directly coupled to the classical environ- 
ment; furthermore, superconductors themselves have a gap such that the 
environment has a low density of low energy modes potentially leading to 
decoherence and the charge fluctuations are screened such that the electric 
coupling to the environment is negligible. The magnetic coupling to the en- 
vironment cannot be completely avoided as the phase changing around the 
frustrated loop is coupled to an induced current, however, the latter can be 
made small if the Josephson junctions are weak and the loop inductance is 
small. In particular, we shall show that the interaction of these currents 
with the environment leads to a small decoherence even if the environment 
includes potentially dangerous normal wires. Finally, we note that the cou- 
pling with the frustrating magnetic field present in the simplest design cau 
be eliminated in designs involving d-wave super conductor sQ or vr— junctions. El 
The first problem to be addressed in the development of a SCPQ is 
the implementation of a convenient probe which tests whether the device 
produces coherent Rabi oscillations. It is obviously difficult to test a sys- 
tem which is decoupled from the environment and it is even more difficult to 
check that it remains in the phase coherent state without disturbing it. Here 
we suggest an experiment to probe the state of such qubits and to determine 
their decoherence time. The idea is to measure the Andreev conductance 



Andreev Spectroscopy for Superconducting Phase Qubits 

between the ground and a normal wire in the geometry of a fork with the two 
prongs connecting to two separate islands of the qubit (Fig. 1). If the qubit 
phase is fixed (classical regime), the low-T conductance exhibits a contribu- 
tion from processes in which electrons from the normal metal are reflected as 
holes from the qubit boundary, diffuse to the bulk superconductor boundary 
and are reflected from it as electrons. Such a phase coherent electron diffu- 
sion in the normal metal prongs of the "fork" leads to periodic conductance 



oscillations (with period $o) as a 



unction of the magnetic flux penetrating 



through the Sf region of the fork.Q'Q This is due to the fact that the mag- 
netic field controls the superconducting phase difference between the two NS 
contacts of the fork [(1) and (2) in Fig. 1] and thus influences the electron 
interference picture. A similar mechanism controls the experiment discussed 
in the quantum regime when the phase of the qubit fluctuates. Qualitatively, 
at small bias voltage the quantum phase acquired by the reflected electrons 
fluctuates and the resulting contribution to the conductivity averages out. 
At larger voltages the phase fluctuations become slow compared with the 
electron tunneling time and the contribution to the conductivity is restored. 
Below we describe the specific device parameters and derive quantitative 
formulas describing this effect. A similar proposal — to probe the Joseph- 
son splitting of levels in a superconducting single Cooper-pair transistor by 
measuring the Andreev conductance — was put forward in Ref. |9[ 



2. SUPERCONDUCTING PHASE QUBIT AND 
ENVIRONMENT 

Consider the SQUID loop with four junctions shown in Fig. 1. Two 
degenerate states naturally appear in such a loop if the flux $g of the external 
field through the qubit loop is exactly $o/2 = hc/Ae. Indeed, in a gauge with 
Ax = ^, Ay = Bx the classical minima of the Josephson energy correspond to 
the phase drops or vr around the loop; in these states the phase of the island 
G (see Fig. 1) is (j) = ib7r/2 (relative to the phase at point 0); below we refer 
to these states as | j) and 1 1). In order to reduce the parasitic coupling to the 
environment, the inductance C of the loop should be sufficiently small so that 
dc/c < W^^^o, Ic is the critical current. In the absence of any charging 
energy in the junctions the system prepared in one of these states will stay 
put forever; quantum effects are due to the small but finite charging energy 
Ec = e^/2C determined by the capacitances C of the junctions, which 
should be smaller than the Josephson energy Ej = hlc/2e » Ec- The 
tunnelling frequency between the two classically degenerate ground states 
is estimated as Qq ~ ^o exp(— a-\//i/cC/e^), where hcoo = ^JSEqEj is the 
Josephson plasma frequency of the contacts; the coefficient a ~ 1.6 was 
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found by a numerical evaluation of the saddle-point trajectory in the space 
of the three phase variables 4'2, 4>3 = 4', and (/>4. We assume values of Ic in the 
range 10~^ A, and capacitances of order of a few fF, resulting in loo/2tt ~ 100 
GHz and a tunneling frequency r2o/27r in the range of few GHz. Our choice 
of parameters is limited by the following constraints: (a) the system should 
be sufficiently small to fulfill the condition CIc/c < 10^'^, but not too small, 
since (b) we need Eq <C Ej] finally, (c) Oq should not be too small, in order 
to keep the decoherence time tdc much longer than ^q . 

Once tunnelling is taken into account, the true eigenstates become | 0) = 
[I T) + I i)]/V2 and | 1) = [| T) ~ I i)]/v^ and are separated by the energy 
gap hi^Q (here and below we shall use the "spin" notation to describe the 
classically degenerate states). The deviation of the magnetic flux $q = BSg 
through the loop from the point ^g = $o/2 removes the degeneracy of the 
states I I) and | |). Thus the intrinsic Hamiltonian of our qubit written in 
the basis | T)) I i) is ^ = hxCTx + hzCTz, where hz = {Ic/V^c)[^q — <l>o/2] and 
hx = hilf). All operations on the qubit can be performed if one is able to vary 
the effective fields hz and hx- The variation of hz can be achieved by changing 
the flux $g in the loop, while the variation of hx can be implemented by two 
means: the variation of the gate potential applied to the island G changes 
this field smoothly (c.f., Ref. [To| ), while switching a capacitor Cext ~ IOC in 
parallel blocks the tunneling abruptly. 

The coupling of a two-level quantum system to the environment is 
usually described by the Caldeira-Leggett "spin-boson" model where the 
"spin" (i.e., tsychJpjel system, TLS) is linearly coupled to an ensemble 
of oscillators .ta'tSliJ Actually, our TLS originates from the dynamics of 
a continuous phase variable (p with a potential energy U{(j)) strongly fa- 
voring values (j) around ib(/>o/2. Then the term that describes the cou- 
pling to the environment can be written, upon integration over environ- 
mental modes, in the form of a non-local imaginary-time action S'diss = 
^ J dtdt'K{t — t')[(j){t) — 0(t')]^. While in a superconducting system we can 
hope the low energy spectrum of the environment to be gapped providing 
favorable conditions for a long coherence time, here, we analyze the (worst 
case) situation where the coupling of low lying modes survives to produce a 
finite ohmic dissipation with K{t) = ri/2'Kt^ , where ?] is the generalized "fric- 
tion coefficient" in the equation of motion r/0 = —dU{(j))/d(j). The strength 
of the environmental coupling is determined by the dimensionless parame- 
ter a = r/0Q/27r^; consider the symmetric system with /i^ = 0: for a > 1 
the tunnelling between the states | t), and | [) is suppressed and the sys- 
tem possesses two classically degenerate ground states. For intermediate 
coupling 1/2 < a < 1 tunnelling is finite but incoherent: at time scales 
t > il^^ = h/hx (with Q the effective tunneling frequency depending on 
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a) the TLS is described by the classical master equation for the occupation 
probabilities P^^^ of the two states with a transition rate T{a) vanishing at 
a ^ 1. At hz = this equation takes the form m{t) = —Tm and the "mag- 
netization" m{t) = {az){t) = Pi{t) — -P|(i) monotonically decays to zero. 
Here, we are interested in the case of very weak coupling ql ^ 1, where 
the TLS dynamics is weakly perturbed by the environment ,Ej i.e., the de- 
coherence time tdc = r~^ shall be much longer than the period of coherent 
oscillations 27ril~^. In this case, the magnetization exhibits damped oscil- 
lations m(t) = exp(— Ft) cos(Ot) with a danming rate T and a frequency 
determined by the parameter a < 1/2 viat3'Ej 

n = {nr/TTa)sm[TTa/{l-a)], T = n/Q = ntan[TTa/2{l - a)], (1) 

with the renormalized tunnelling frequency O^ = r2o(i^o/'-^c)" ; here, 
u>c is an upper cutoff for the frequencies of the environmental modes. A 
measurement of m(t), particularly its quality factor Q, then allows us to 
test the degree of quantum coherence in the device. 

The direct determination of m{t) in a time-domain experiment like 
the one reported in Ref.ll£l (where a superconducting "single Cooper-pair" 
device was studied) is a demanding task. As an alternative one may 
attempt to fmd the correlation function D(t) = {az{t)az{0)) through a 
measurement^ of the rf absorption; indeed, the latter provides the imag- 
inary part x"{^) of the response function x{^) = 2 J^ dtD"'{t) eyip{iujt), 
where D°'{t) = {[az{0),az{t)])/2i is the antisymmetric part of the spin-spin 
correlation function. Standard relations then allow us to extract the correla- 
tor D(t) from the measured absorption x"i^)'- the symmetrized correlation 
function ^^(t) = {{azi0),crz{t)})/2 = coth{hu;/2T)x"iuj) relates to the cor- 
relator via D{lu) = 2D%uj)/[l + exp(-u;/r)] (note that D{u; < 0) = at 
T = 0). Below we discuss yet another experiment giving us access to D{t): 
we will see that the nonlinear Andreev conductance dlj^/dV through the 
qubit's island G can be expressed via the Fourier-transform D{(jj) of its cor- 
relator. Scanning magnetic field and voltage in a transport measurement 
then appears to give a much easier and direct access to the operation of the 
qubit than the more usual real-time and r/-techniques. 

Being not aware of exact results available for xl"-^)— ,or D^i^) = 
{{az{0),crz{t)})/2, we make use of an approximate analysisEJ to arrive at 
an idea of the relevant dependencies in these quantities: for Q ^ 1 the 
absorption characteristic is close to a combination of Lorentzians: x"i^) = 
F/[F^ + {lo — r^)^] — F/[F^ + (cj + rj)^], implying a symmetrized correlator of 
the form D-%uj) = 4\Lo\Qr/[{r^ + Q"^ + uj'^f - Aoj'^Q?]. In the time domain, 
D'^{t)\T=o ~ cos(Qt) exp(-Ft) + A/[i]:Q{Q,t)\ where the last term is due to 
the non-analyticity of D''{uj) at a; = 0. 
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3. ANDREEV CONDUCTANCE 

We first concentrate on the subgap conductance through a NIS contact 
(7^1 in Fig. 1) with a low normal-state conductance at = Gte^ /h <^ e'^/h and 
discuss the effect of fluctuations in the phase (j){t) on the island G. Consider 
an ideal SCPQ operating at the degeneracy point {hz = 0) with the ground 
and excited states |0) and |1) separated in energy by HQ. The operator 
C^ = exp(— i(/)) creating Cooper pairs changes the relative sign of the states 
1 1) and II), thus C^\0) = |1) and C' has no matrix element in the ground 
state. Therefore |0) is a state with a well-defined number of Cooper-pairs 
modulo 2, implying that the usual Andreev process transmitting a single 
Cooper pair through the interface is forbidden at low voltage V < hi}/2e, 
while at voltages above the gap hQ/2 the Andreev conductance will be large. 
Higher-order processes with simultaneous tunnelling of two Cooper pairs are 
allowed, but their rate is of the order of Gf and thus is negligible at Gj <C 1. 
Overall, we then expect an Andreev conductance dlj^/dV similar to the one 
of a usual NIS junction but with the superconducting gap A replaced by the 
two-level separation fiVl. On the other hand, for a fixed phase on the island 
we can expect to observe the usual Andreev current across a NIS junction. 

The above consideration defines the simplest example of a situation 
where the quantum fluctuations of the superconducting phase suppress the 
Andreev conductance of the SIN contact. General expressions describing 



these effects have been derived in Ref. 21 (c.f.. Appendix B). Here, we discuss 



the result for the nonlinear subgap current /a(^) for the case where the 
normal metal is relatively clean (with a dimensionless conductance g ^ 1) 
such that all corrections of relative order \/g can be neglected, 

X [tanh[(£; - E' /2 + eV)/2T] - tanh[(£; + E' /2 - eV)/2T]] . 

Here, G{E) is the real part of the Cooperon amplitude in the normal metal 
(see below) and D^{E) = J dt exp{iEt/h)D^{t) is the Fourier-transformed 
autocorrelation function D^{t) = (^(t)^"'"(0)) of the superconducting order 
parameter ^{t) = exp[i(/)(f)] on the island G. On long time scales t ^ 
r2~^ the oscillatory modes near the minima of the Josephson potential are 
irrelevant; the values of (j) in these minima differ by tt, thus D^{t) is equivalent 
to the (sought for) TLS correlation function D{t) discussed in the previous 
section. 

The real part of the Cooperon amplitude is given by 

C{E) = -^R y — —, r , (3) 



i^V ^ Dq^-2iE/h + T^^' 
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where the sum goes over the Cooperon eigenmodes g„, V denotes the volume, 
v the density of states (for a single projection of spin), and D the diffusion 
constant of the normal metal. Assuming a wire geometry with w, L ^ 
{uiin{LE = y/hUjAE, L^ = y^DT^),w), and d < w denoting its length. 



width and thickness, Eq. (y) reduces to the expression C{E) ss pLe/[1 + 
(?i/4£'T(^)^)]^'^, where the resistance p = {2vDwd)~^ per unit length of the 
wire is measured in units of (/i/e'^)/Length. The term oc r~^ describes the 
electron dephasingoresent at all T > due to electron-electron interactions 
and external noise.E3 For a short wire, C{E) = pL; in the end we arrive at the 
simple result C{E) « pm.m{LE, L,n^). An expression for Ia{V) analogous 
to Eqs. (^) has been derived earlieicHI for a similar physical problem. 

Let us analyze the result @ in the limit T — > which applies to the 
regime T <^ (eV, hQ). At the sametime we keep the electron dephasing rate 
T^^ nonzero, since experimentally" an apparent saturation at T < 0.1 — IK 
has been found. Differentiating Eq. (^) by ^ at T = we find 



dh 
dV 



271%^ 



2eV 



dED^{E)C{eV -E/2) 



(4) 



(values E <Q are excluded from the integral since D^{E < 0) = at T = 0). 
Consider first of all the degenerate case with /i^ = in the absence of 
any dephasing within the qubit, i.e., E = and D^{E) = 2'k5{E/Ti — Q). 
The evaluation of the integral in (0) is trivial and we obtain the result 



dl^ 
dV 



0, eV < hn/2, 

{e^/n)GlC{eV -hVL/2), hVL/2<eV, 



(5) 



in agreement with the above qualitative arguments. In the following we re- 
strict the discussion to long wires L » y^JT/Ti; then the conductivity above 
threshold shows a square-root singularity which is smeared by a nonzero de- 
phasing rate r~^. The nonzero conductance below threshold is due either to 
a finite dissipation rate E > or a finite temperature. At zero temperature, 
a finite E produces the sub-threshold {eV < hQ,/2) conductance 



diA 
dV 




eV <. nn/2, 



6 = nn-2eV <. nn/2, 



where we have assumed r^ — > oo and have made use of the form of D(j){E) 
discussed at the end of Sec. 2. The numerical coefficient a depends on the 
shape of the correlator D^[E) and is of order unity; for the Lorentzian shape 
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discussed above we find a = 32/37r. A finite phase breaking time r^ in the 
normal wire changes this result at low voltages eV < fi/r^: dlj^/dV ~ 
{e^ /fi)GlpL^p[liT{eV)'^ /{h^Y]. Finally, a finite temperature also produces 
a finite sub-threshold conductance which, however, can be well separated 
(at T <^ C^eV, hO, — 2eV)) from the above behavior due to its exponential 
oc exp[(2ey — hQ)/T] rather than algebraic dependence oc {■vKTeV/fiQ,)^ 
and oc hT/{hQ — 2eV) at low and high voltages. Note, that these results for 
the suppression in the conductance dlp^/dV due to phase fluctuations refer 
to terms of order Gl and hence our analysis is valid only for low Gt ■ 

The above discussion clarifies the appearance of a threshold behav- 
ior in the degenerate case h^ = 0, where the differential conductance 
dlj^/dV sharply drops below the gap hQ/2 and vanishes at F, T — > 0. 
For hz > hx = hQ the correct ground state is close to a state with a 
well-defined phase (p 3.nd the conductance attains its semiclassical value 
dl'^/dV = {e^ /2K)G1 pyJfiD / eV . As a consequence, a sharp minimum (at 
$q = $0/2 corresponding to h^ = 0) is expected in /a(^) when changing the 
magnetic flux ^q through the qubit. This qualitative conclusion survives the 
presence of a finite electron decoherence rate r~^ in the normal-metal, how- 
ever, the precise form of dlj^/dV is modified. A short electron decoherence 
time Tip < h/T is frequently observed in thin wires at lowest temperatures 



as relevant for the experiments on qubits (c.f., Ref. 23). These observations 
are poorly understood theoretically and thus it is more practical to mea- 
sure T^ independently in the same experimental setup. Such a measurement 
can be carried out via an analysis of the magnetic field dependent Andreev 
conductance in the "fork" geometry as we are now going to discuss. 

Consider the interference experiment for the device shown in Fig. 1. 
The first stage of the proposed experiment is to measure the amplitude 
I^ in the oscillations of the Andreev current /a = /a + -^a + -^a ^ a-s a 
function of magnetic field B with the aim to extract the electron decoherence 
time T^p in the normal metal (the superscripts '*^ refer to the contributions 
from the contacts i = 1,2, see Fig. 1). The total phase determining the 
interference current I^ sums up to i;^a = i;^ + 27rBSf/^o. Away from the 
degeneracy field Bn = (n + 1/2) ^o/Sg, the (semiclassical) phase (j) in our 
device is determined by the minimization of the Josephson energy and we 
obtain (p = Tr{{BSg/^o + 1/2} — 1/2), where {x} is the fractional part of 
X. Choosing a geometry with Sg = Sj, we find for BSg/^o £ [n,n + 1/2] 
an interference current Ia^' = 2ji2 cos (f>A = 2ji2 cos{37rBSq/^o), whereas for 
BSg/^o G [n 1^1/2'^ + 1] the sign changes: /I'''' = -2j^2Cos{^7rBSq/^o)- 
The amplitude3 of the interference component J12 depends on the precise 
geometry of the "fork" and on the electron dephasing length L,p which should 
be longer than the distance L12 between the contacts #1,2 (otherwise J12 will 
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1/2 BSg/% 1 



Fig. 2. Schematic view of the Andreev current through the "fork" enclosing 
the qubit versus magnetic flux for the special case Sq = Sj. The fork serves 
two purposes: the interference pattern gives access to r^ through measuring 
I^ and allows for the determination of the exact point of degeneracy. 



be exponentially small). Detailed calculations and measurements of ji2 have 
been presented by the Saclay grouppHSIHSI they found a discrepancy by a 
factor 1.8 between the best fit (with r^ as fitting parameter) to their data 
and the calculated expression for ji2 (cf-, Ref. 20, p. 190; this discrepancy 
may be due to the use of a simplified description of the electron decoherence 
using a factor exp(— t/r,^) in the real time expression for the CooperonB'Ej). 
Second, we test for the coherent time evolution of the qubit: In the geometry 
with Sq = Sf, the semiclassical current amplitude I]^ (B) is periodic with 
a period ^o/Sq, showing upward cusps at B = Bn, see Fig. 2, dotted line. 
These cusps are due to the penetration of single flux quanta $0 into the 
qubit loop at B = Bn- At these fields, the interference contribution I}^ 



r(2) 



vanishes and the semiclassical Andreev current is given by the sum 1)^+11 
(in the generic case Sq 7^ 5/ the component I}^ does not vanish at B = Bn 
and the flux penetration leads to a jump in the measured Ia). In addition, at 
B = Bn the quantum fluctuations in the phase (j) become large and suppress 
the Andreev current li to the "active" island G of the SCPQ, see Eq. (Q). 
Hence, close to the points B = Bn a narrow dip is superimposed upon the 
above-mentioned cusp in the Ia{B) dependence, see Fig. 2, solid line. The 
relative width AB/Bq ~ eQ/Ic of this dip is determined by the relation 
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hz = Ic{^q — *^o/2)/\/2c < hx = hO, and provides the coherence gap of the 
qubit. In addition, we can measure the decoherence rate T via an analysis 
of the conductance diA/dV = o"a(1^) at the degeneracy points i? = S„. In 
the hmit T <C eV and eV — > the entire conductance is determined by the 
current through the contact #2, o"a'(V^) = {e'^/h)G^C{eV). The variation 
of the full differential conductance ^^{V) = (Ta'(^) +o"a'(^) at low voltages 
eV < hQ then can be used to extract the value of T. Note that C{E) 
for the fork differs from the form used in the expression for the subgap 



conductance, c.f., Ref. 19. The detailed analysis of (Ta(^) for the different 
relations between fi, F, and r^ will appear in a future publication. 

Finally, any method measuring the SCPQ coherence time is useful only 
if the presence of the measurement circuit does not, by itself, lead to a de- 
coherence rate comparable to the "intrinsic" one. In our setup there are two 
sources of additional phase (p decoherence due to the presence of the normal- 
metal wire nearby: (i) the direct coupling of the qubit to the environment 
via the non-zero Andreev conductance, and (ii) the inductive coupling be- 
tween ac supercurrents in the qubit loop and normal currents in the Andreev 
fork. We consider them one by one: (i) Comparing the dissipative action 
for the Andreev conductance (c.f., Ref. ^ for example) with the general 
relation for Sdiss and the definition of the parameter a = r](pQ/2iTh one finds 
the contribution Oa = {'JTh/8e'^)af, where af = {e"^ /fi)G^pin.m{L,L^) is 
the semiclassical Andreev conductance via the contact #1. Thus, a classi- 
cal subgap resistance of the SN contacts in the MJ7 range will produce a 
value a A < 10~^. (ii) We calculate the friction coefficient r/ind via the en- 
ergy dissipation rate Wind = iTmA{d4> / dt)"^ . The main contribution to Wind is 
due to the ac magnetic field applied to the conducting wire (with thickness 
d < 50 nm, width w ~ 100 nm and a few microns length L), which gives 
Wind ~ ^{wLd)\" {uj)B'^, where ^"{oj) ~ ojw'^ /c^ is the imaginary part of the 
magnetic susceptibility of a thin-film wire in a magnetic field B^ in the low- 
frequency limit (for frequencies uj ~ lO^'' and a conductivity a ~ 10~^ ficm 
of the wire the effective skin depth is much longer than the wire's width 
w ~ 100 nm). The magnetic field outside the SCPQ loop (of size Lq ~ l//m) 
is less than ^q/L"^, thus (ujE^)"^ < {uj(f)^)'^{h/e)'^{^q/^o)'^L~'^. Combin- 
ing all factors, one finds aind < {w'^ L/ pLq){^q/^of . For a SCPQ with 



^q/^o < 10 ^ and w/Lq ~ 0.1 this results in a negligible value Oind < 10 ^. 



'q/ ^u ^ -L" ""'•■^ uji -^q ~ 0.1 this results 

In conclusion, we have derived expressions for the subgap liy) char- 
acteristics, Eqs. (|2|) and (Q), which can be used for the determination of 
the intrinsic correlation function D^f, = (exp[i^(t)] exp[— i(/)(0)]) of a super- 
conducting phase qubit (SCPQ), once the electron dephasing time r^ that 
enters the Cooperon, Eq. (^), is known. As the latter is difficult to control 
theoretically at the ultralow temperatures of the experiment we propose an 
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interference experiment probing the SCPQ with an "Andreev fork" , allowing 
to determine r^ within the same setup and then use its value to extract tdc- 
In the absence of decoherence in the SCPQ, F = 0, the current /a(^) is zero 
(at T = 0) below the threshold eV = hi}/2, irrespective of the value of r^. 
Thus the SCPQ decoherence time tdc = P^^ can be measured even if it is 
longer than Tip. In practice, the most important limitation of this technique 
seems to be due to the non- vanishing temperature of experiment which pro- 
duces a nonzero sub-threshold conductance even at P = 0. However, in the 
range T <^ fiQ and at voltages such that T <^ {2eV, hQ — 2eV) the behavior 
of dlp^/dV is exponentially close to the zero-T one; therefore even values of 
P <C T/fi can be extracted from the measured dlj^/dV curves. 
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for helpful discussions and the 'Zentrum fiir theoretische Studien' at ETH- 
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Science via the programm "Physics of Quantum Computing" . 

REFERENCES 

1. P. Shor, Los Alamitos Conference Proceedings, 1994. 

2. L.K. Grover, Phys. Rev. Lett. 79, 325 (1997). 

3. A. Ekert and R. Jozsa, Rev. Mod. Phys. 68, 733 (1996). 

4. L.B. loffe, et al, Nature 398, 679 (1999). 

5. J.E. Mooij, et al, Science 285, 1036 (1999); |cond-mat/990828S . 



6. G. Blatter, V.B. Geshkenbein, and L.B. loffe, to he published. 

7. Yu.V. Nazarov and F.W.J. Hekking, Phys. Rev. Lett. 71, 1625 (1993). 

8. H. Pothier et al, Phys. Rev. Lett. 73, 2488 (1994). 

9. F.W.J. Hekking, L.I. Glazman, and G. Schon, Phys. Rev. B 51, 15312 (1995). 

10. D.A. Ivanov and M.V. Feigel'man, JETP 87, 349 (1998). 

11. A.V. Veretennikov et ai, Proceedings of LT-22. 

12. S. Chakravarty and A.J. Leggett, Phys. Rev. Lett. 52, 5 (1984). 

13. A.J. Leggett et ai, Rev. Mod. Phys. 59, 1 (1987). 

14. U. Weiss, Quantum Dissipative Systems, World Scientific (1998). 

15. F. Lesage and H. Saleur, Phys. Rev. Lett. 80, 4370 (1998). 

16. Y. Nakamura, Yu.A. Pashkin, and J.S. Tsai, Nature 398, 786 (1999). 

17. Y. Nakamura, CD. Chen, and J.S. Tsai, Phys. Rev. Lett. 79, 2328 (1997). 

18. K. Volker, Phys. Rev. B 58, 1862 (1998). 

19. H. Pothier, S. Gueron, D. Esteve, and M.H. Devoret, Physica B 203, 226 (1994). 

20. S. Gueron, PhD thesis, Saclay (1997). 

21. M.V. Feigel'man, A.I. Larkin, and M.A. Skvortsov, |cond-mat/9907358| . 



22. A. Huck, F.W.J. Hekking, and B. Kra mer, Europhys. Lett. 41, 201 (1998). 

23. M.E. Gershenzon, cond-mat/9908099|, Localization-99 Proceedings. 



24. M.V. Feigel'man and A.I. Larkin, Chem. Phys. 235, 107 (1998). 



